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ABSTRACT 

A gravitational observatory such as LISA will detect coalescing pairs of mas- 
sive black holes, accurately measure their luminosity distance and help identify 
a host galaxy or an electromagnetic counterpart. If dark energy is a manifesta- 
tion of modified gravity on large scales, gravitational waves from cosmologically- 
distant spacetime sirens are direct probes of this new physics. For example, a 
gravitational Hubble diagram based on black hole pair luminosity distances and 
host galaxy redshifts could reveal a large distance extra-dimensional leakage of 
gravity. Various additional signatures may be expected in a gravitational signal 
propagated over cosmological scales. 



1. Introduction 

Evidence for an accelerating expansion of the U niverse is getting s tronger. A mea- 
surable dimming of dist ant type la super novae (e.g., lAstier et al.l l2006l : iRiess et al.l 12007 



Wood-Vasev et al.lboOTh. CMB data (e.g.. ISpergel et al.l 120071 ) and additional cosmological 



probes (e.g., IWrightl 120071 ) indicate that a repulsive "dark energy" comprises 73% of the 



Universe's content. Characterizing the properties of dark energy and uncovering its physical 
nature are two of the most important goals of modern cosmology. Current observational 
strategies include tests of the p ossibility that dark e nergy arises from a failure of general 
relativity on cosmological scales (lAlbrecht et al.ll2006l ). 



Essentially all astronomical measurements are performed via electromagnetic waves. 
The availability of accurate gravitational wave measurements, within the next decade or so, 
will thus be a significant development. In particular, since the propagation of photons and 
gravitons could differ at a fundamental level, we argue here that gravitational waves emitted 
by cosmologically-distant "space-time sirens," such as coalescing pairs of massive black holes, 
may be used as valuable alternative probes of dark energy physics. 
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Black holes with masses ^, IQ^Mp, are present at the ce nter of numerous nearby galaxies 



(e.g. iKormendy fc Richstond Il995l : iMagorrian et al.l Il998l ). As such galaxies collide over 
cosmic times, their central black holes coalesce, releasing > 10^^ ergs of binding energy in 
the form of gravitational waves (hereafter GWs). To measure the GWs emitted by these 
cosmologically-distant space-time sirens, ESA and NASA will build the Laser Interferometer 
Space Antenna, LISA0. 



2. Gravitational Measurements 



GWs emitted by black hole binaries have the unusual property of providing a direct 
measure of the luminosity distance, D^, to the black ho les, without extrinsic calibration. 
Owing to the highly coherent nature of GW emission (jSchula Il986l ). the amplitude (or 
strain), /i+x, frequency, /, and frequency derivative, /, of the leading order (quadrupolar) 
GW inspiral signal scale as 
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where +x represents the two transverse GW polarizations, Mc = {mim2Y^^ / {mi +m2)^^^ 
is the black hole pair "chirp" mass and z its redshift. Provided the GW source can be 
reasonably well localized on the sky, an extended observation of the chirping signal leads 
to precise measurements of /i+x, /, / and thus Dl, independently. LISA's orbital con- 
figuration allo ws for a "tr i angulation" of GW sources on the sky to within a solid angle 
5fi ~ 1 deg^ (ICutlerlll998l : IVecchidl2004l ). thus providing very accurate dis t ance measure- 
ments, with SD jDr < 1% at ^ < 2 typically JOutlei] [l998l : iHughesI 12002 : Ivecchiol bood : 
Lang fc Hugh^booeh . 



Recently, the possibility of identifying the individual host galaxy in whic h a pair of merg- 
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Kocsis et al. 


2006, 


2007) 



and identifications out to ^ 5-7 may even be possible in some cases (IKocsis et al.l 120071 ). 
A unique host galaxy identification can be achieved through coordinated observations with 
traditional telescopes, either to survey the LISA-triangulated area for unusual galactic prop- 
erties or central activity after the coalescence, or by monitoring in real time the sky area for 
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unusual electromagnetic emission, as the coalescence proceeds. T ypically >, 10 ergs of ki- 
netic energy are delivered to the recoiling black hole remnant (e.g.. lSchnittman fc Buonanno 



20071 . and references therein) and its environment. The disturbed gas surrounding coalescing 



black holes could thus power bright electromagnetic emission during and after coalescence 
(lArmitage fc Natarajanll2002l : iMilosavljevic &: Phinneyl l2005l : iDotti et al.ll2006l ). permitting 
the coincident identification of a unique host galaxy. 



3. Gravitational Hubble Diagram 



A consequence of successfully identifying the host galaxies of coalescing black hole pairs 
is the possibility to draw a gravitational Hubble diagram, i.e. one that relates the gravita- 
tional luminosity distances, Dl, of space-time sirens to the electromagnetic redshifts, z, of 
their host galaxies. 

The solid line in Fig. [T] shows a model fit to a Hubble diagram from a recent type la 
supernova sample with a dark energy equation of state P/ p = w = —1.05 ±0.3 (at 3a, shown 
by the shaded region). Luminosity distances are expressed in terms of the usual distance 
modulus, u = 5 logm f- PfV 1 Mpc) + 25. The individual supernova data points of the gold 
sample of iRiess et al.l (120071 ) are also shown. A fiat universe with Hq = 72 km s~^ Mpc~^ 
and a matter density flm = 0.27 is assumed. 

The general interest of a gravitational Hubble diagram can be illustrated by considering 
alternative explanations for the apparent dimming of distant type la supernovae. If this 
dimming were caused by dust atten uation, photon-axi o n conversions or an i ntrinsic evolution 
of th e supernova population (e.g., ICsaki et al.l 120021 : iMirizzi et al.l l2006l : iRobaina fc Cepa 
20071 ). rather than stretched distances, this would become apparent in comparisons between 
gravitational and electromagnetic Hubble diagrams. Indeed, in each of these alternative 
dimming scenarios, photon- and graviton-based Hubble diagrams would be fundamentally 
discrepant since self-calibrated gravitational distance measurements are not susceptible to 
any significant bias from absorption, scattering, reddening, or axion- conversion. 

In practice, however, the value of such a comparison is limited by line-of-sight matter 
inhomogeneities, which generate "weak lensing" uncertaint i es on any photon or graviton 



D]^ measurement (IHolz fc Hughes! l2005l : iKocsis et al.l l2006l : iDalal et al. 



20061 ). While the 



lensing effect can be averaged out over the many random lines-of-sights available with the 
supernova data (i.e., data points in Fig. [1]), it may not be the case for coalesci ng black hole 



pairs if LISA merger event rates are modest (e.g., a few tens per y ear at z < 5; iMenou et al. 



Wvithe & Loeb 
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individual measurements amount to dis tance uncertainties ranging from SDl/Dl — 1% at 



0.5 to 5Di/Di ~ 10% at 2; = 5 (e.g. lKocsis et al.ll2006l ). The corresponding Sex distance 



modulus uncertainties exceed the Sa confidence contours on the current dark energy model 
fit (shaded region in Fig. [T]) at z > 0.7, which makes comparisons between photon- and 
graviton-based Hubble diagrams imprecise even at moderate redshifts. The extent to which 
LISA events can be used to perform meaningful comparisons with photon-based Hubble 
diagrams will thus depend strongly on the actual distribution of massive black hole merger 
events with redshifts and the corresponding efficiency of host galaxy identifications. 



4. Modified Gravity 



The greatest prospect for dark energy science with gravitational waves may lie in explor- 
ing new physics on cosmological scales. The possibility that accelerated expansion results 
from a failure of general relativity has fueled much theoretical work on large scale modifi- 
cations of gravity over the past few years. Since building a satisfactory theory of modified 
relativistic gravity is a formidable task, any insight that can be gained from direct observa- 
tional constraints on the linearized gravitational wave regime cannot be overlooked. 

One may expect gravity modifications to contain a new length scale, let us call it Re-, 
beyond which gravity deviates from general relativity. In order to explain the observed ac- 
celerated expansion of the Universe, this scale is expected to be of the order of the current 
Hubble radius H^^ . Modified gravity must also pass standard tests of general relativity on 
scales much shorter than Re-, e.g. in the solar system and in the strong field regi me of binary 



pulsa rs. An existence proof of modifications of this type is given by DGP gravity (IDvali et al. 



2OOOI ). a braneworld model with an infinite (possibly many) e xtra dimension. In this model. 



which also leads to an accelerated expansion of the Universe (jPeffayet II2OOII : iDeffayet et al. 
2002a), gravity is intrinsically higher dimensional^. To leading order, the gravitational po- 
tential has the standard 1/D behavior at distances D smaller than i?c; while it behaves 
5-dimensionally at larger distances (1/D^). Moreover, the stronger the gravitational field, 
the closer the theory is to general relativity. As a consequence the model passes the PPN 
tests in the solar system, \ yith possible deviations emerging in upco ming generations of so- 
lar system measurements (ILue fc StarkmanI l2003l : iDvali et al.l 120031 ) . The G W emission of 
classical astrophysical sources is also expected to closely match that of general relativity. It 
is only at very large distances and low curvature that gravity is beginning to "leak" in the 



Gorbunov et a 



and 



( 2006h . lCharmousis et all |2006h . lDeffavet et al.l (|20oi ) . IPvalil |2006h llzumi et al.l|2007h 



Gregory et al.l (|2007( ) discuss the stability of the DGP self-accelerating phase. 
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extra dimension. 

A previously unexplored consequence of extra-dimensional leakage is that cosmologically- 
distant GW sources would appear dimmer than they truly are, from the loss of GW energy 
flux to the bulk. Inspired by this idea, we investigate here the consequences of this possible 
modification of gravity, borrowing from the DGP model the notion that strong field grav- 
ity (and hence GW emission) asymptotes to general relativity, while deviations appear at 
very large, typically cosmological distances, in the weak field regime. We will not deal here 
specifically with GWs in the DGP model, which is the subject of a separate study, rather, 
we illustrate more generally how a gravitational observatory such as LISA may reveal cos- 
mological deviations in the weak field graviton propagator. 

In the presence of large distance leakage (say at distances much larger than Rc), flux 
conservation over a source-centered hypersphere requires that the GW amplitude scales with 
distance D from the source as 

where dim is the total number of space-time dimensions accessible to gravity modes. Thus, 
for dim > 5, it deviates from the usual h^^{D) cx 1/D scaling. The scaling in Eq. 
i s consistent with explicit GW calculations in spac etimes with compact extra-dimensions 



dimensions open up only at large distances (iGregory et al.l 120001 : iDvali et al.ll2000l . l200ll ) 



(ICardoso et al.ll2003l : iBarvinsky fc SolodukhinI 120031) and also applies to models where ex tra 



The top three lines in Fig.[T]show gravitational Hubble diagrams, i.e. the expected locus 
of Dl and z measurements from black hole pairs and host galaxies, in three simple scenarios 
with dim = 5 and leakage beyond a scale Rc obeying Eq. ([3]). To allow arbitrary possibilities 
for the transition at the cross-over scale, we adopt /i_|_x(-Dl) oc {Dl[1 + {Dl/ Rc)"^'^Y^'^)~^, 
where n determines the transition steepness. In the three cases shown, cross-over scales of a 
few Hubble distances, Rc = 1-4 Dh ~ 2-9 Gpc were used and values = 1 or 10 ("steep") 
for the transition steepness were adopted. Error bars show the magnitude of la weak lensing 
uncertainties on Dl measurements due to line-of-sight matter inhomogeneities. 

The electromagnetic Hubble diagrams of each scenario is assumed to mimic the cosmol- 
ogy of a w = —1 dark energy model (sohd line), as expected if these scenarios produce an 
accelerated background cosmological FLRW space-time which affects the geodesic motion 
of photons in the standard way. As a result of the different propagation of GWs and large 
distance leakage of GW energy flux, however, measured GW amplitudes, /i+x, are reduced 
and gravitational values are overestimated (Eqs. [l]-[2]) for sources beyond the cross-over 
scale in these scenarios. The corresponding graviton Hubble diagrams would thus deviate 
from the electromagnetic version (solid line) since all photon-based measurements are stan- 
dard, by construction. The discrepancy between the graviton-measured and photon-inferred 



- 6 - 



Dl values could be of order unity even at moderate redshifts. Black hole merger events and 
associated host galaxies may thus reveal the leakage of gravity with cross-over scales ^ a, 
few Hubble distances. 

Diagnostics based exclusively on electromagnetic measuremen ts cannot reveal gravi- 



tational leakage directly, although t hey may do so indirectly (e.g., iLue fc StarkmanI 12004 



Wang et al. 


2007; 


Zhane et al. 


2007) 



identify with only a GW signal. As is well known, all the mass combinations that determine 
the GW signal, e.g. in post-Newtonian expansions of the inspiral s ignal, are redshifted by 



isp 

3 SI £ 



a factor 1 + z and thus degenerate with redshift (e.g., lHughesll2002l ). Provided that general 
relativistic black holes and GW emission are accurately recovered on sub-cosmological scales 
in the modified gravity scenarios considered, leakage on cosmological scales would lead to 
a hidden bias on the redshift (overestimated) and black hole masses (underestimated), as 
inferred from adopting a background cosmology that ignores leakage altogether. Only if the 
redshift of the host galaxy were to become available would the discrepancy with the grav- 
itational measurement become apparent, as illustrated in Fig. [H Paradoxically, this could 
limit our abilities to identify the host galaxy of a merging black hole pair, by invalidating 
search strategies that rely on a (biased) valu e of the gravitatio nal luminosity distance to 



define redshift cuts on potential host galaxies (iKocsis et al.ll2006l ). 



We note that evidence of GW leakage could also emerge in different contexts, e.g. from 
l ow redshift short gamma-ray bursts detected by a network of ground-based GW detectors 
( iDalal et al.ll2006l ). or from an unexpectedly low amplitude of the present-day GW back- 
ground relative to its CMB- normalized value (A. Buonanno, priv. communication). 

Despite the strong GW leakage suggested by the above energetic arguments, one should 
be very cautious in evaluating the significance of the deviations shown in Fig. [H Indeed, 
Eq. ([3]) says very little about the cross-over transition physics and other details of GW 
propagation on cosmological scales. In particular, even though th e scaling given in Eq. ([3]) 
does apply to the DGP model, the "infrared transparency" effect (IDvali et al.ll200ll ) can be 
shown to result in a considerable increase of the distance at which GW leakage is manifested 
i n DGP, reaching sc ales much beyond Rc for sources with frequencies relevant to LISA 
( iDeffayet et al.l 120071 ). Therefore, the DGP model does not necessarily produce significant 
deviations between gravitational and electromagnetic luminosity distances and it would, in 
fact, be difficult to distinguish from general relativity on the basis of the test illustrated in 
Fig. [TJ On the other hand, since it is presently unclear whether infrared transparency and 
its consequences for GW propagation are generic to all cases of higher- dimensional gravity 
with large distance leakage (in particular when Lorentz symmetry is broken). Fig. [1] remains 
useful in providing a good measure of potential GW leakage. 
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5. Additional Signatures 



Large distance leakage is only one of several possible modified gravity signatures in 
the GW signal from cosmologically-distant spacetime sirens. We discuss a few additional 
possibilities here, voluntarily adopting a simple phenomenological approach. 

A first class of signatures resides in the GW polarization signal. In many modified grav- 
ity scenarios, additional polariz ations exist beyond the two transverse quadrupolar (-|-x^ 



modes of general relativity (e.g., IWill 
theor ies (including the f(R) va riety; 



20011 ) and vector-tensor theories ( iBekenstein 



2006h. This i s the case, for instance, iri scalar-tensor 



Wands 



1994 



2004 



Maggiore fc Nicolis 



Jacobson fc Mattingly 



2000 : iNakao et al. 



20041 ). Braneworld 



gravity scenarios are no exception since the Kaluza-Klein graviton extra polarizations can 
also in principle be radiated by various types of sources. The no-hair theorem is traditionally 
invoked to forbid any extra polarization modes in the signal from two coalescing black holes. 
However, it is well known that black hole unicity th eorems fail to apply in the usual sense 
to higher dimensional theories (lEmparan et al.ll2002l ). Moreover, in DGP gravity the gravi- 
ton extra polarization is expected to show up at large d istances even for (possibly static) 
spherically symmetric space-times (jPeffayet et al.ll2002bl ). Black holes can also be hairy in 
theories of massive gravity (where the graviton carries extra polarizations with respect to 



those of a massless graviton), in particular when Lorentz invariance is broken (IBlas et al. 
20061 . 120071 : iDubovsky et al.l 120071 . and references therein). 



A second class of signatures is related to the GW signal propagation velocity which, in 
modi fied gravity scenarios, c an differ from the speed of light. Propagation can be subluminal 



[e.g.. IDubovsky et al.ll2005l ) or superluminal (e.g.. I Jacobson fc Mattinglyl |2004| ) . The pos- 
sibility to time a GW signal propagated over cosmological distances, relative to the signal 
from a prompt electromagnetic counterpart causally associated with the black hole merger, 
may thus offer additional diagnostics of large-scale modified gravity. In addition, signatures 
may emerge because the number of spacetime dimensions available to gravity is odd. In this 
case, even a massless mode ha s a Green's function which does not vanish inside its light cone 
see e.g., ICardoso et aP booA Consequences for GW pro pagation in braneworld gravity 
have been explored in part (IBarvinsky &: Solodukhirul2003l ). but not fully in the case when 
gravity is modified at large distances. Thus, the oddness of spacetime on cosmological scales, 
like in DGP gravity, could add to the set of modified gravity signatures. 

Finally, a third class of signatures relates to the phase of the GW signal, which could 
deviate from general relativistic expectations once propagated over cosmological distances. 

In summary, GWs from cosmologically-distant spacetime sirens may be valuable alterna- 
tive probes of modified gravity since various signatures may exist in a GW signal propagated 
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over cosmological distances. For some of these signatures to become apparent, the identifi- 
cation of a host galaxy or an electromagnetic counterpart to the spacetime siren is required. 
Conversely, the absence of such signatures would constrain gravity modifications and provide 
a consistency check on other methods employed to discover the nature of dark energy. 



This work has made use of the advanced cosmology calculator (jWrightl 120061 ). It is a 
pleasure to thank A. Buonanno, G. Esposito-Farese, G. Gabadadze, Z. Haiman, D. Helfand 
and B. Kocsis for useful inputs and discussions. 
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Redshift z 

Fig. 1. — Electromagnetic Hubble diagram of a standard dark energy model (solid line + 
shaded region) and hypothetical gravitational Hubble diagrams for three modified gravity 
scenarios with biased gravitational luminosity distances from large distance gravitational 
leakage, beyond a cross-over scale Rc (top three lines). 



